Universal Features of the Time Evolution of Evanescent Modes in a Left-Handed 

Perfect Lens. 
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The time evolution of evanescent modes in Pendry's perfect lens proposal for ideally lossless and 
homogeneous, left-handed materials is analyzed. We show that time development of sub-wavelength 
resolution exhibits universal features, independent of model details. This is due to the unavoidable 
near-degeneracy of surface electromagnetic modes in the deep sub-wavelength region. By means of 
a mechanical analog, it is shown that an intrinsic time scale (missed in stationary studies) has to be 
associated with any desired lateral resolution. A time-dependent cut-off length emerges, removing 
the problem of divergences claimed to invalidate Pendry's proposal. 
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Long ago, Veselago pointed out that very unusual 
properties, such as negative refraction, would be exhib- 
ited by materials with negative refraction index. Re- 
cently, Pendry has claimed that those so called left- 
handed (LH) materials with e(cjo) = A*(k>o) = — 1 ; can act 
as perfect lenses with, ideally, arbitrary sub-wavelength 
resolution. In addition to its genuine conceptual impor- 
tance, this proposal has attracted much attention due to 
the practical realization of man-made materials expected 
to be left-handed Q , where negative refraction has been 
claimed to be observed Q . The challenge to conventional 
ideas conveyed in Pendry's work has fueled a heated de- 
bate that has contributed to sharpen the issue, if not to 
settle it HUSH 

The criticism expressed by Garcia and N-Vesperinas |7j 
(see also Refs. 0,0) seems of particular importance, for 
it would imply that some sort of fundamental violation 
of physical laws is unavoidable in Pendry's perfect lens. 
The idea is that the necessary amplification of evanescent 
modes would turn a square-integrable incoming wave into 
a non normalizable signal, somethingdeemed unaccept- 
able by the authors of Refs. [1013. Pendry has 
replied that losses, unavoidable in the real world, would 
provide a natural cut-off, forbidding the amplification of 
large wavevectors (though degrading the perfect resolu- 
tion). Haldanc 10] has put the blame of divergences on 
the role played by surface modes (polaritons) local- 
ized at the interfaces of the LH material: amplification 
of evanescent modes is a gift of those modes |2j. For a 
homogeneous material, those surface modes exist and be- 
come dispersionless for increasing wavevectors, being the 
unphysical absence of an intrinsic cut-off length which 
causes the pathologies [ljj- Any realization of a LH ma- 
terial (think of the composite nature of proposed LH ma- 
terials Q| or photonic crystals [l^) must have a natural 
length below which the homogeneous description fails, 
providing the necessary cure to divergences at the price 
of lower resolution. 

Although Pendry's and Haldane's escapes from diver- 
gences are certainly safe, they seem to suggest that the 



textbook idealization of lossless and homogeneous me- 
dia, that works so well otherwise [lflj . is fundamentally 
flawed when applied to a material that happens to satisfy 
e(u>o) = /Lt(cJo) = — 1; a t some frequency ujq. Apparently, 
only after the inclusion on the real world constraints 
(losses and/or small-distance structure) could the pro- 
posal be made physically acceptable. I find this state of 
affairs very unsatisfactory, and it is the purpose of this 
paper to show that, even within the self-imposed ideal- 
izations of a lossless (for u) — luq) and purely homoge- 
neous, left-handed material, Pendry's perfect lens pro- 
posal is correct. 

To this end, we will consider the time evolution (l4Hl5| 
of sub- wavelength features. Time development requires 
the study of the frequency dispersion of a particular 
model (or material) , apparently preventing us from draw- 
ing general conclusions. However, it will be shown that 
the structure of the relevant magnitudes in the immedi- 
ate vicinity of the target frequency is dominated by the 
surface polariton modes [llj, . These modes show uni- 
versal features for any LH material at long wavevectors, 
therefore allowing us to extract general conclusions. A 
key point will be the identification of a time scale for 
any intended length resolution. This time scale, linked 
to the near-degeneracy of surface modes and missed in 
stationary studies, will provide the clue for overcoming 
the problem of divergences. 

Consider a transverse current oscillating at frequency 
u>, located in a source plane at x = xq (see Fig. HJ. The 
Fourier component (k) of the current leads to an electric 
field given by (S polarization, for simplicity): 



E=[0,E k (x,t),0] 
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for evanescent 
E k (x,t) ~ e -p\x-xa\ e -i<»t } w here p = Vk 2 - uj 2 , in 
units with light velocity c = 1 (the connection between 
current and field need not concern us here [lfij)- In 
the presence of a slab of LH material between x = ±o, 
the field acquires the usual reflected and transmitted 
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FIG. 1: Geometry of the problem: slab of LH material be- 
tween planes x = ±a, with source plane at xo, and image 
plane at xi. The curve of Ek(x) illustrates perfect restora- 
tion of evanescent fields when xq = —2a and x\ = 2a. 

components (only the vacuum part explicit and for 
arbitrary normalization of the incoming wave): 

E k (x, t) = [ e -P\*-*°\ + r k {uj) e P( x + x ^} e~ luJt , x < -a 
E k (x, t) = t k (uj) e -P( x - x °) e~ luJt , x > a 

(2) 

Pendry [2| has shown that, for a frequency ojq such that 
e(u>o) = (J,(u>o) = — 1, the reflection and transmission co- 
efficients are: 

r k (aj ) = 0, t k (uj ) = exp(4p a) (3) 

with 

po = aA 2 -c 2 (4) 

This leads to a perfect restoration at the focal plane 
X\ = 2a of a source field at xq = —2a, as shown in Fig. 
□ 

The exponential amplification inside the LH material, 
necessary for image restoration, is at the root of the di- 
vergences pointed out before 0, Q . Notice that a square 
integrable field with Fourier components £ k at the source 
plane xq — —2a, will emerge at the plane x = a with 
norm: 

N= f d\ l \E(x = a,r ll )\ 2 ocJ2\^e Poa \ 2 , (5) 

k 

This norm certainly diverges for large wavevectors, unless 
unjustified constraints are put on the source field. 

So far we have described the stationary solution, but 
our aim is the time development of such a situation. 
Therefore, the frequency behavior of r(u>) and t(u>) is 
needed. In the vicinity of u> = u/q, and in the deep sub- 
wavelength region, the reflection and transmission coef- 
ficients exhibit a highly singular behavior controlled by 
the appearance of simple poles. This singular behavior 
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FIG. 2: Left panel: Real (continuous line) and imaginary 
(dashed line, artificially enlarging delta functions) parts of 
the transmission coefficient ifc(w) in the vicinity of Uo. Right 
panel: dispersion relation of surface modes u)k,±, illustrating 
near-degeneracy for large wavevectors. Insets show the spatial 
form of both modes in the weak-coupling regime. 

can be described by: 

(6) 

and is depicted in Fig. [21 (left panel). One can see that, 
for instance, t k {uj) is basically the sum of four simple 
poles located at frequencies ±u>k t ±. Notice that the tar- 
get frequency is sandwiched between two such poles 0J k ,±, 
corresponding to the surface modes (right panel of Fig. 
|2J). These are the even and odd combinations of the fun- 
damental mode of an isolated interface, that would take 
place exactly at u> = wq. The frequencies u>k,± are solu- 
tions of the equations [llj: [tanh( y o m a)] ±1 = —p(uj)p/p m , 
with p m = ^Jk 2 — e(uj) p(lo)ui 2 . This gives: 

w2 i± =wg±»£ (7) 

with an exponentially small coupling 10] for large 
wavevectors 

4^8Cu; 2 e- 2 ^ (8) 

where (Cc^o) -1 = 2p' + (e' + ^)(fc 2 /w 2 - with 
fi' = dn(u>o) / dcu and e' = de(u>o)/du>. It is important to 
realize that the expressions of Eq. EJ though only valid 
in the neighborhood of ojq, do contain the singular struc- 
ture associated with the surface modes. Therefore, the 
relevant dynamics of these modes is expected to be well 
reproduced. Ruppin has studied these modes thor- 
oughly, and their importance for the amplifying process 
has been noticed before [3 fiol . 

We can now study, for instance, the time evolu- 
tion of the transmitted field at the right interface 
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Etrans(x = a,t), in terms of the incident field at the left 
interface Ei nc {x = — a,t). This is enough to understand 
the amplification issue. We will assume a pure sinusoidal 
wave at the frequency ujq, but with a well defined time 
origin, chosen to coincide with the signal arrival at the 
left interface: 

E mc (x = -a, t) = 6{t){E Q e^ at + c.c), (9) 

9(t) being Heaviside's unit step function. 

The evaluation of E tra ns{x = a,t) is now trivial but, 
given the unusual and controversial nature of this sub- 
ject, we choose to change the language in the hope of 
bringing the results to a far more familiar situation. Let 
us say that we have two identical (left and right) oscil- 
lators, (xi,x r ), with natural frequency u> , and a weak 
coupling ijI that splits the degeneracy w| ± = uj 2 ± r] 2 , . 
Now we force the left oscillator with the external force 
f(t) and watch the dynamics of the right oscillator. The 
equations of motion are: 

xi + jii + u)%x t + rfex r = f(t) , 1Q 
x r + 7± r + u 2 x r + rilxi — 

where a damping 7 has been added for later convenience 
but, in accordance with our idealization, is supposed to 
be 7 = 0, for the moment. Upon identifying the ex- 
ternal force, left, and right oscillators with the incident, 
reflected (E re f), and transmitted fields in the following 
manner: 

/(*) «— » e 2po ° E inc {x = -a,t) 
x\ < — > E re f(x = -a,t) (11) 
x r < > Et rans (x +&,£), 

this simple problem is entirely equivalent to our origi- 
nal one. The interpretation of Eq. ^] is direct: the 
incoming wave plays the role of an external force hitting 
the left interface, and exciting a reflected (left oscillator) 
and a transmitted (right oscillator) wave. This mapping 
allows us to understand the physics of the original situ- 
ation in simpler terms. For instance, forcing the system 
(left oscillator) with frequency w , the stationary solution 
tells us that, surprisingly, only the right oscillator moves. 
This corresponds (through Eq. Ill|) to the absence of a 
reflected wave. 

The time development of the field corresponding to the 
incoming perturbation of Eq. [5] is now easily obtained, 
with the following result: 

E trans (x = +a, t) = A(t) e 2poa 9{t)E e-^ ot + c.c. (12) 

where a characteristic time-dependent modulation am- 
plitude A{t) appears: 

= 1 - (e- |Aw +* + e-**"-*) / 2 (13) 

where Aoj± = cjfc,± — — ±Awfc/2 for |Awfc| << luq, 
with the splitting of polariton modes Alo^ given by: 

Au k = uj k , + - ~8Cuj e- 2paa (14) 



in the deep sub-wavelength region. This beating of 
modes, familiar from the mechanical analog, is sketched 
in Fig. |31 Notice that the total response is the superpo- 
sition of the stationary solution (o>o), which is Pendry's 
solution, with the normal modes (w±). The latter are 
unavoidably excited and do not decay in time owing to 
the absence of losses around uj . 

The relevance of this behavior for the problem of di- 
vergences should be clear by now. For a fixed lapse of 
time t (large in units of the bare period: uj^t >> 1, but 
otherwise arbitrary) , short wavelengths corresponding to 
splittings smaller than Auik ~t~~ x , have barely begun to 
emerge at the right interface. Therefore, we can identify 
a time-dependent crossover wavevector k e ^(t), satisfy- 
ing the condition: 

t Au) k * f} = 1 (15) 
such that, for k >> k e ^ , then 

A{t) - (Au; k t) 2 - (CW) 2 e~ 4poa (16) 
The norm of Eq. is now replaced by 

AAoc (Cu;ot) 4 J2\^ 3p ° a \ 2 , (17) 

fe 

for fc >> k e ff . Therefore, amplification has been re- 
placed by decay, solving the problem of divergences. 
Of course, progressively shorter details take longer and 
longer time (and more energy from the current source) 
to develop 0, 0] , but nothing pathological affects the 
physics of the system. I emphasize again that this picture 
is unavoidable in the deep sub- wavelength limit, irrespec- 
tive of the model details of e(w) and fi(to) ■ It amounts 
to recognizing that, for a given wavevector fc, the sys- 
tem's dynamics has a characteristic time scale given by 
the (inverse) of the frequency splitting between the two 
surface polaritons corresponding to that wavevector (Eq. 
1141 ) . The relevance of this time scale would remain hid- 
den if only stationary solutions were analyzed. It is a 
curious though rewarding fact that, in a system that 
has been idealized to be homogeneous down to arbitrary 
small scales (absence of intrinsic length scale) , time takes 
up the responsibility of providing the necessary cut-off 
length in order to prevent divergent catastrophes. 

It could be argued that the solution provided here to 
the problem of divergences is an artifact of the ideal 
lossless nature of the problem. Real systems would 
show losses that will eventually kill any transient regime, 
perhaps leaving again the divergent stationary solution. 
This is not the case, as can be seen by the explicit inclu- 
sion of a finite life-time (r = 7 _1 ) in Eq. A stationary 
solution is indeed established, but with amplitude A sta t 
greatly reduced from the lossless limit A sta t{oo) at large 
wavevectors: 

Astatjr) _ 1 / 18 x 

A stat (oo) l + (rAw fe )- 2 
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t / (2 jr / Ato k ) 



FIG. 3: Time evolution of the modulation amplitude (in 
units of the stationary, nearly lossless, limit) for three values 
of surface life-time corresponding to zero (r Auik = oo), weak 
(r Au>h = 10), and strong (r Auik = 0.5) damping. 

Notice that, for small values of tAw/j, the stationary am- 
plitude is what would have been expected if the lossless 
evolution were suddenly stopped at a time of the order of 
t. Therefore, we see exactly the same physics as in the 
lossless case, but now the role of the observation lapse 
of time is taken by the polariton life-time. A resolution 
wavenumber (k e ff) can be defined again (compare Eq. 

T AWfee,/ = 1 (19) 

such that, for wavenumbers larger than this cut-off, ex- 
ponential decay rather than amplification is observed at 
the exit of the slab. The norm of field at the exit is now 
as in Eq. 1171 with the exchange t <-> r. 

The transient behavior for the two regimes 
t >> (Aujk)^ 1 and r < (Aw^) -1 is illustrated in 
Fig. The picture is simple: when the life-time of 
surface polaritons is much longer that the characteristic 
time scale (inverse splitting of modes) for the wavevector 
k, the stationary limit approaches the lossless case. This 
means that surface modes have had enough time to 
build the final response before dying away. On the other 
hand, if the surface modes do not live long enough, no 
amplification is possible. 

In the deep sub- wavelength limit, Eq. El allows us to 
provide an explicit expression for the minimum life-time 
(Tmin) required to get a resolution l res = 2ir/k e ff , with 
radiation of wavelength Aq = 2tt/loq: 



Train ~ exp(47ra\/ l T es - Aq 2 ) (20) 

o G luq v 

Similar results have been obtained before for the resolu- 
tion , further reinforcing the correctness of our restric- 
tion to the polariton dynamics. Notice the characteristic 
exponential dependence on resolution. Although mainly 



concerned with matters of principle in this paper, this de- 
manding result clearly shows that the polariton life-time 
may well be a major limiting factor for a practical realiza- 
tion of perfect lenses. I believe this vulnerability of 
the ideal situation, a fingerprint of the near-degeneracy 
of surface modes, is at the root of similar results found 
before [HOH. 

In spite of potential difficulties in the road to a practi- 
cal left-handed amplifier, the essential physics described 
in the paper may have been observed in a different sys- 
tem. The role of coherence between surface modes (and 
the required time scale associated with it) in the expla- 
nation of the extraordinary transmission through hole 
arrays [19| is strikingly similar to our treatment. This re- 
inforces the view that both problems (left-handed am- 
plification and extraordinary transmission through hole 
arrays) are probably different manifestations of the same 
physical behavior of surface modes. 
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